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Abstract. Systems of classical continuous particles in the grand canonical ensem- 
ble interacting through purely attractive, yet stable, interactions are defined. By a 
lattice approximation, FKG ferromagnetic inequalities are shown to hold for such par- 
ticle systems. Using these inequalities, a construction of the infinite volume measures 
is given by a monotonicity and upper bound argument. Invar iance under Euclidean 
transformations is proven for the infinite volume measures. The construction works for 
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Q\ • arbitrary activity and temperature and for integrable long range interactions. Also, 



inhomogeneous systems of particles with different "charge" can be treated. 
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1 Introduction 

In the thermodynamics of continuous, classical particles interaction is usually 
defined through pair potentials. This however does not allow to model purely 
attractive interaction, as some kind of repulsive core has to be included into the 
pair potential in order to render it thermodynamically stable. Similar consider- 
ations hold for potentials including n-point interactions with some upper bound 
on n. 

At the same time, the repulsive core up to now is an obstacle to gain better 
control over Gibbs measures outside the low-density — high-temperature (LD- 
HT) regime. For this regime Gibbs measures have been constructed rigorously 
by various kinds of cluster expansions [T3], see also [3| for a LD-HT cluster ex- 
pansion of a class of non-pair potentials. With the exception of [121 EH] using 
Pigorov-Sinai cluster expansion for low temperatures, the results on Gibbs mea- 
sures outside of this regime mostly use Banach-Alaoglu like compactness argu- 
ments in combination with tightness to establish existence, see e.g. jH |Hl ITU] . 
Both methods however do not give complete control over the number of limit 
measures even for a specified set of boundary conditions. 
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New applications of random measures and point processes in "soft matter 
physics", sociology and biology have furthermore increased the need for models 
that are not based on the rigorous assumption that the interaction is defined 
through a pair potential with repulsive core, cf. e.g. [T71 l2Uj . 

Constructions of Gibbs measures for arbitrary parameters have been obtained 
using GKS ferromagnetic inequalities for (discrete) spin systems or (continuous) 
two-dimensional Euclidean quantum field theories [18j . It is the aim of this paper, 
which is based on (21 E] , to apply such Euclidean quantum field theory methods 
in connection with FKG inequalities |5] to continuous particle systems. FKG 
inequalities for continuous particle systems have been considered before, see [7] 
for a general discussion with applications to Widom-Rowlinson and continuous 
random cluster type models. Even though in that reference there is no discussion 
of the thermodynamic limit, a number of estimates used for the construction of 
this limit in Section 4 below can, in the special case of only one particle type, 
also be formulated in terms of the stochastic ordering criterion, which is the main 
result of [7j, see 0. In this work, a different method which also extends to more 
than one type of particles and is based on comparison through interpolation, is 
used. 

The paper is organized as follows: In Section 2 we define a new class of purely 
attractive, yet stable, potentials for systems of continuous, classical particles in 
the grand canonical ensemble. These are Kac-like potentials given by concave, 
and linearly bounded energy densities in the static field generated by charged 
particles. Such potentials model attractive forces between particles of the same 
type and repulsive forces between particles of different type. They could e.g. 
be useful to describe the attractive part of van der Waals forces of water and 
hydrocarbon molecules. The linear bound, which is necessary for stability reasons, 
then implies that the difference in energy of e.g. a water molecule completely 
immersed in hydrocarbon molecules and a water molecule completely immersed in 
water molecules is given by a chemical potential, which is reasonable. Compared 
with non-stable purely attractive pair potentials the saturation of attractiveness 
between particles of the same type for the potentials used here none the less can 
be seen as a weak attractive behavior, which explains the title of this article. 
On the physical side, these models for the special case of two particle types with 
opposite charge (or spin) have some similarity with ferrofluid models [51 [TU|. fTT| . 
The latter models however use repulsive forces between particles of all types 
to obtain stability and are mathematically rather different from what is being 
discussed here. 

In Section 3 ferromagnetic FKG correlation inequalities are proven through 
an approximation with lattice gas or, in other words, with spin systems. These 
inequalities are then used to construct the themodynamic limit (TD limit) first on 
the level of expectations of monotonically increasing function and then on the level 
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of the measures describing the ensemble, see Section 4. The construction works for 
arbitrary temperature and density (activity), in particular no low density - high 
temperature conditions are needed. The infinite volume measures are uniquely 
determined and turn out to be invariant under Euclidean transformations. The 
space dimension plays no role and also long-range interactions can be treated as 
long as they are integrable. 

2 Weakly attractive interaction 

Let us begin with the heuristics of weakly attractive interaction. For simplicity, 
in the heuristic discussion we only consider one type of (identical) particles in 
the Euclidean space M. d with d G N the space dimension. We assume that each 
such particle carries a positive unit charge. A charge s G R in the point y G M. d 
in the point x G M. d gives rise to a static field sG(x — y), where G is a symmetric 
{G(x) = G(—x)) non-negative function in L 1 (M d , dx). The static field <fi generated 
by a finite number n of particles is given by <p = G * r] where rj = J^jLi ^% where 
Uj G M. d , j = 1, ... ,n give the positions of the particles and 5 y stands for the 
Dirac measure in y. Let v : R — > R be the energy density of the field (p. For 
technical reasons we assume that v (0) = 0, v G C 2 (R) and v has bounded first 
derivative, |f'(0)| < b V0 G R and some constant b > 0. The potential energy of 
the particle distribution rj is given by integration over the energy density of the 
associated static field 



Let now v be concave, i.e. v" < 0, cf. Fig.l, and y G R be the position of 
a further particle giving rise to the field A<p(x) = G(x — y). For the heuristic 
discussion we assume that G is monotonically decreasing in the distance. Then, 
if y is far away from the support of rj, the field $ ext ' = G * rj can be neglected 
in the "range" of G around y. The extra particle in y at a point q in this range 
of G around y then leads to a change of the energy density by Aej = v(A<f)(q)). 
Suppose now the points y and q are translated towards the support of n where 
the "exterior" field $ ext ' can no longer be neglected. The difference in the energy 
density at the translated point q' is then given by Ae2 = v (<3> ext '(g') + A0(g)) — 
v ($ ext '(g')). As v is concave and $ ext ' > 0, we get Aei > Ae2, see Fig. 1. Hence, 
the particle in y gains energy if it approaches the support of r\. This explains 
the attractive character of the potential introduced in (|TJ for concave v. The 
potentials of this type are called weakly attractive potentials, because, due to the 
linear asymptotics of v , there is a saturation of the attractiveness for the single 
particle in regions which are homogenuously filled by a very strong "exterior" 
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field $ ext - ; see again Fig. 1. If we allow the particles in yi,...,y n to carry 
charges s\, . . . , s n G R, hence r) = YTj=i s j fiy 3 , = G * r]. It is then clear from 
the above discussion that weakly attractive potentials, as defined by (0) with v 
concave, lead to attractive interaction between particles of the same charge and 
to repulsive forces between particles of different charges. 

In a thermodynamic ensemble of particles with weakly attractive interaction, 
groupings of particles with the same sign of the charge are energetically preferred 
w.r.t. mixtures of particles with different signs of the charge. This leads to posi- 
tive correlation between the event of finding a particle of given charge in a certain 
region and the event of finding another particle nearby with the same sign of the 
charge. In technical terms, this should be expressed by a "ferromagnetic" correla- 
tion inequality. On the other hand, the saturation of attractiveness in regions of 
strong fields of either sign implies that the number of particles in a given volume 
can always be controlled (also in the TD limit) by the number of particles of a 
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free gas with higher activity. These two ingredients - ferromagnetic correlation 
inequalities, implying monotonicity of the expectation values for monotonically 
increasing function in the TD limit, and upper bounds of the particle number 
due to comparison with a free gas of higher activity - is essentially all what is 
needed to copy Nelson's classical proof of the TD limit of P ((^-theories |18j . 

Having clarified the heuristic basis, we now carry on with the technical formu- 
lation following [21 El- Let z > be the activity and let r be a probability measure 
on K with compact support, supp r C [-C, C] for some C > 0, such that r{0} = 0. 
We define a purely Poisson Levy characteristic ip{t) = z f R [e lst — 1] dr(s), t E M, 
and we set for / G <S(R d ), where <S(M d ) is the space of real Schwartz test functions 



Then, Cq : S(WL d ) — > C is a positive definite, normalized functional which is also 
continuous w.r.t. the nuclear Schwarz topology on iS(M d ). Let 5'(R d ) be the 
space of tempered distributions over M. d (i.e. the topological dual space of iS(M d )) 
and let B = B(S' (IR d )) be the Borel sigma algebra (the sigma algebra generated 
by open sets) over iS'(IR d ). By Minlos' theorem QjJ there exists a unique 
probability measure fi on the measurable space (S'{R d ),B) such that 



The measure /Iq has a natural interpretation as the Gibbs measure of a non- 
interacting system of classical, continuous particles in the grand canonical en- 
semble where each particle carries a random charge distributed according to r. 
In fact, for A C M d with finite Lebesgue volume, |A| < oo, let be a Pois- 
son random variable (number of particles in A) with intensity z|A| and {Yj}j G fq, 
{Sj}j e N two families of i.i.d. random variables (also mutually independent and 
independent from N^) where Yj is uniformly distributed in A (position of the jth 
particle in A) and Sj (charge of the jth particle in A) has distribution according 
to r, then the restriction of the coordinate process (.,/) to A, which maps test 
functions with support in A to random variables on (S' (R d ) , B , fi ) is equivalent 
in law with the random process S(A) 3 f — > X,j=i Sj f(Xi) = (Vo,a, f) where 



For the measure /x this together with a short estimate [3] implies that /io(r) = 
1 where T = T c C 5'(R d ) is the set of all signed measures with only finitely many 
support points in any compact subset of M d s.t. for rj G T the following holds: 
\v{y}\ < C \/y E R d and J Rd (l + \x\ 2 )-^ d+e ^ 2 d\r]\(x) < oo Ve > 0. For rj = 
E^esupp^ s v 6 v, \V\ = J2y & u PPV \ s y\ 6 v T is in fact a measurable set in (S'(R d ), B), 
cf. 0. In the probability space (iS'(lR d ), B, /io) one can replace 5'(R d ) with T and 
B with B r = B H T, the trace sigma algebra of B on F. 




C (f) =e^ d ' W) 



(2) 




(3) 
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These support properties of /aq imply the following: Let the kernel function 
G be as above. One can show that for 77 G T, G * 77 G Ll oc (Mr , dx) [Sj- For a 
infra-red cut-off function with compact support g G Co(lR d ), # > 0, one can thus 
define the following interaction 

U g {n)= / v(0)gdx, (f> = G*r). (4) 

Here t> fulfills the conditions given above, in particular v is concave. That we 
defined the infra-red cut-off in this way instead of using finite volume measures 
/xo,a and the interaction (JTJ) is only a technicality, which is going to allow us to 
keep as close to Nelson's proof of the TD limit for P{4>)2 as possible. 
From j2] we now get: 

Proposition 2.1 (i) U g : T — > K is B r measurable; 
(11) U g G LP(T, fio) Vp > I; 
(in) e- u » G L p (r,/i ) Vp> 1. 

Sketch of the proof: The crucial point in the proof of Proposition 12.11 is the 
very simple estimate 

\U 9 ( V )\<b [ G*gd\rj\ (5) 

JR d 

and the observation that 

J e b{lvl ' G * 9) dLi (r]) < 00 (6) 

where the latter inequality follows from calculating the Laplace transform of the 
measure llq o tp^ 1 with f+(r]) = |r/|, cf. (j2^jl below. For the details see [3J. ■ 

By Proposition 12.11 (iii) the following interacting Gibbs measures 11 g with 
infra-red cut-off g are well-defined on (r, B r ): 

e -u g (v) 
J r e u b dfio 

The aim of the paper is now to remove the infra-red cut-off g, i.e. for inverse 
temperature (3 > we want to investigate the limit g /* (3 for g G Co(IR d ) for the 
measures \x g and show their convergence to a limit measure Lip. Before we can 
resume with this task in Section 4, ferromagnetic FKG inequalities have to be 
established through lattice approximation, which is done in the following section. 
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3 Lattice approximation and FKG inequality 

First let us introduce the FKG (for Fortuin, Kasteleyn and Ginibre) inequalities 
and state the main result of this section: 

Let 0(T) be a specific set of functions - also called observables - i.e. mea- 
surable functions F : T —>■ R. Namely, for F G 0(T) there exists n G N, a con- 
tinuous function H : R n — > R exponentially bounded and fast falling functions 
h u ...,h n G Cf. f .(R d ) such that F(r/) = #((77, /ii), . . . (77, fc*)). The fast falling 
functions are defined as C f . f (R d ) = [h G C(R d ) : sup xgRd |(1 + \x\ 2 ) N h(x)\ < 00 
ViV G N} and exponentially bounded means that there exist k, K > s.t. 
-£/~(xi, . . . , x n ) < Ke K (\ Xl \ + '" + \ Xn ^ Wx G R n . From the definition of T it is clear 
that (r],h) is well-defined for r\ G T, /i G Cf.f(R d ). By the exponential bounded- 
ness of the functions H, \F(r})\ < Ke K ^' h \ h = \ht\ + ■■■ + \h n \ G C f . f .(K d ) and 
by an estimate like © it follows that 0(T) C LP{T,fx ) Vp > 1. Proposition O 
(iii) also implies 0(T) C L p (r,/i 9 ), p > 1. 

An observable F G C(r) is called monotonically increasing if F can be rep- 
resented by a function H : M n — > R and h\,...,h n G C f . f .(R d ) such that hi > 0, 
I = 1, . . . , n, and H(x\, . . . , x n ) is monotonically increasing in each argument 
Xi,...,x n . By 0(r) we denote the collection of all monotonically increasing 
functions F G 0(T). By definition, a probability measure /1 on (r, B r ) fulfills the 
FKG inequality, if 

J F X F 2 dn> J F x dyL J F 2 dn VFi, F 2 e 6{T). (8) 

Theorem 3.1 The FKG inequality (0) holds for ri g defined as in Section 2. 

The proof of the FKG inequality for the gas of continuous particles in the 
GCE with weakly attractive interaction requires various steps of approximation 
preserving the FKG inequality. How the approximations are done in detail to a 
large extend is a matter of convenience. Here we do the following steps: 

• Approximate the continuous system with kernel G G L 1 (R d , dx) as specified 
in the previous section by a continuous system with finite range interactions 
G G C (R d ), G > (Proposition O; 

• Approximate the continuous system with the regular and finite-range kernel 
G by infinite lattice systems (Proposition 13. 3|) : 

• Approximate infinite lattice systems by finite lattice systems (Proposition 



Verify the FKG inequality for the finite lattice system (Proposition 13. 5|) . 
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The technical details of the different approximation always follow the same 
scheme, namely a proof of convergence in L p (r,/Jo) through pointwise conver- 
gence and a L P (T, /j )-upper bound through an estimate as in the proof of Propo- 
sition To avoid repetitions, we give the detailed proof for the first steps only 
and then indicate where this scheme has to be modified for the remaining steps. 

Let G be as described in the previous section and let G n G C (M d ), G n > 0, 
and G n — > G in L 1 (IR d ,(ia;) as n — > oo such that \G n \ < G for a L 1 (IR d ,(ia;) 
majorant G. Let ?7 ff(n be defined as in (jlj) with G replaced by G n . Then 

Proposition 3.2 (%) C/ gjn — > U g in L p {T,fi ) as n — > oo /orp > 1; 
(^zj e _c/9 ' n — > e -17 ® m L p (r,/io) as n ^ oo /brp > 1. 

Proof. We notify that 

\U 9 (v)-U g , n (v)\ <b [ \G-G n \*gd\ V \, V eT (9) 

G n -y G in L 1 (M. d , dx) and # G C (K d ) imply |G - G n \ * g -> pointwise in M d as 
n — » oo. By (JHJ) , L d G * (7 d|7/| < 00 for /i almost all 77 G T. Thus, the right hand 
side of (jHJ) converges to zero by the theorem of dominated convergence applied 
on L 1 (IR d , d\rj\) for //q almost all rj G T as ri — > 00. 

This implies that C/ 9) „ — > U g and e _c/ff '" — > e~ Uf! /io a.s. as n — > 00. The 
assertion (ii) now follows from the theorem of dominated convergence in L P (T, /j,q). 
In fact, e~ Ug ' n ^ < e b ^ ,G * 9 ' > uniformly in n and the right hand side of this estimate 
is in L p (T,fi ) for p > 1, cf. (jHJ) where 6 can be chosen arbitrarily. The assertion 
(i) follows from the fact that the upper bound in (ii) can always be taken as an 
upper bound for U gtn . ■ 

In the remainder of this section we thus assume G G Cq(M. ), G > 0, if not 
defied differently. The next step is the lattice approximation: 

Let A > be the lattice spacing for the lattice XL d C M. d . We consider each 
point of this lattice as center of a box A^, j G AZ, which has side length A. 
The boundaries of the boxes are attached to either of the boxes having these 
boundaries in common, such that the give a disjoint partition of W d . Given 
r] G T, we define rjj = t](Aj). Define rj x = J2je\z d ^ dr lj 1a a with 1a the indicator 
function of the set A C M d .The lattice interaction Ug : V — > R is defined through 

U g K (v)= [ v(<p x )gdx, ^ = G* V X . (10) 

By a simple adaptation of the proof of Proposition 12.11 to the discretized case 
one immediately gets that the statements (i)-(iii) of Proposition 12. II remain true 
if U g is replaced by U x . Furthermore, for F G 0(T) defined through F(j]) = 
H((rj, hi), . . . , (77, h n )) as described above, we obtain F x by replacing rj by r/ A . 
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That this procedure might depend on H and hi,...,h n and not only on F, is not 
an obstacle: For any F G 0(T) we a priori fix one such representation. We get 
the following convergence of the lattice potential and the lattice observables: 

Proposition 3.3 (i) U x -> U g in L p (T,u, ) as\\0forp> 1; 
(it) e~ U 9 — > e~ Ua as A \ in L p {T,fi Q ) for p > 1; 
fmj F A -> F as A \ in L P (T, ri ) forp>l. 

Proof. Let h G Cf. f .(K d ) and rj ET, rj — J2 y& u PPV s v V Then ' 

= ^ \~ d r]j j h(x)dx 
je\i d A j 

»/ A. 

j£\Z d ^GsupprynA^ J 

= V s y \- d f h(x) dx — > (77, /i) as A \ 0. (11) 

^ /a A 

yesupp?7 jA(„) 

Here j A (y) for y G suppry is defined to be j G AZ s.t. y G A A . The sums all 
converge absolutely by definition of V and Cf.f.(R d ). The convergence as A \ 
follows from X~ d J AA h(x) dx — > h(y) and a dominated convergence argument 

j x {y) 

using the rapid decay of h and \s y \ < C. 

We first prove (i) and (ii). Note that <p(x) = G * rj{x) = (rj,G x ), G x {y) = 
G(x — y), and G x G Cf.f.(R d ) and the same holds for (f> x (x) with rj replaced by r/ x . 
By (fTTjl. A (x) — > Vx G R d . Consequently, v((j) x ) — ► v(4>) in L x (R d , g dx) as 

|w(0 A (x))| < feC||C||ootl(A(gr, G, A) n supp?]) VxGsuppfif (12) 

with b the linear bound of v , C from the definition of T and A(g, G, A) = {y G R d : 
inf xg ( SU p Pg+SU p pG ) |y — a; | < A} and the term on the right hand side is monotonically 
falling in A. §A denotes the number of elements in a set A. Thus, U x -> U g 

pointwisely in T. A L P (T, /i ) -upper bound for [7" A , A < 1, and e _c/ s is e aNA( - G ' s '^ 
with a = bC\\G\U\g\\ LHRdydx) , N^) = ^AHsuppr/) for A C R d , cf. JEJ. 

Let us now consider (iii). F x — > F pointwisely on T is a consequence of 
(jll|) and the continuity of if : ~R dn — > R. The L p (r, yU ) -upper bound ex- 
ists because of the exponential boundedness of H and the estimate \(r] x ,h)\ < 
Bsup xEMd |/i(x)(l + |x| 2 ) 7V / 2 |(|r/|,(l + |.| 2 )- 7V / 2 ), for N > d, A < 1 and B > large 
enough, where the exponential of the r.h.s. is in L p (r,/io) for p > 1, cf. © and 
0- ■ 

The next step is the approximation of r] x by rj\ = J2je\z d nA ^ d Vj^-A x as 
A C R d compact approaches R d . The notation A R d is used for this limit; the 



10 



H. Gottschalk 



precise meaning of lim A is that for any sequence A n C R d compact such that 
A n C A m for n <m and U rae NA„ = R d the limit exists and all such limits coincide. 
Ug fjjf) = J Rd f(0 A ) gdx, <t>\ = G* rjjl, and F A is defined through a representation 
H ': R dn —> R, hi, . . . , h n G C LL (R d ) by replacing (77, hi) with (7$, h), I = 1, . . . , n. 

Proposition 3.4 (i) U X A -> U x in L p (T,fi ) as A / R d forp> I; 

(11) e~^ A -> e"^ as A / M d m L p (r, /i ) /or p > 1; 
fmj F£ -»• F A as A / M d in L p (T,fi ) }orp > 1. 

Proof. Firstly we notify that the upper L P (T, /j )-bounds given in Proposition 
EH for £/ 9 A e"^ and F x are upper L p (r, /i )-bounds also for U X A , e U ^ A and F A 
uniformly in A C R d (strictly speaking this holds for A < 1 only but related 
bounds can be found for A > 1). It remains to prove pointwise convergence. 

To prove (i) and (ii) we use the finite range of G and we obtain A (x) — > (j) x (x) 
as A y R d where the limit for fixed x G R d is obtained increasing A in the limit 
A / R d finitely many times. Thus, v (0 A ) — > v((f) X ) on R d pointwisely. With the 
same L 1 (R d , gdx)-upper bounds as in Proposition 13.31 this implies U x A — > U x on 
all r. 

It remains to prove (iii). The convergence {r]\,h) (r] x ,h) for 77 G T, h G 
Cf.f.(lR d ) as A /" R d , is a simple consequence of the definition of these spaces 
which permits to get the statement by the theorem dominated convergence in 
L 1 (R d ,dx). This implies the pointwise convergence _F A — > F x on T as A /* R d . 
■ 

Until now we have reduced the system to a finite lattice system where the state 
space can be identified with R m with m = m(A, A) = jjL, L = L(A, A) = A n \Z d . 
Before proceeding to the FKG inequality, we want to make this more explicit. 
In fact, for F G 0(T) with fixed representation ^(77) = H((r), hi), . . . , (77, h n )). 
Thus, keeping A C R d compact and A > fixed, 

F X ( V ) = H\ Y, *- d V*(lA>,hi),..., £ A-% A (l A ,,/i n ) J 
\jeAn\z d jeknxi d J 

= H( V f\j€h). (13) 

We notify that for F G 0(T) one can choose H G <D(R m ) where the latter space 
is defined as the space of exponentially bounded, continuous functions from R m 
to R monotonically increasing in each argument Xj G R, j = 1, . . . , m. For 
different iable functions H this follows by differentiation w.r.t. n x using the chain 
rule and A _o! (1 a a, hi) > if hi > 0. For only continuous H this follows by 
approximation with differentiable functions. 
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We note that the family of random variables {r]j}jeh is i.i.d distributed. Let 
p be this distribution on R, i.e. p is the unique probability measure on (1 



with Fourier transform f R e tst dp(s) = J r e V ' A i dp (r]) = e xd ^\ cf. (J2J) and 
note that tjj(0) = 0. We can now define the lattice measure p^ on (R m ,£>(R m )) 
through 

e -u^{vj\jeh) 

d^j&Jlj) = 7 -uHn UCDj ®m ( -dp® m ^i^m) , (14) 



with U^fa) = U* A (r]) for rjj = 7]}, j e L, i.e. U^irjjlj 6 L) = f Rd vfa) gdx, 
These definitions imply 

H^ Jf^ ( 1B ) 

for any H : R m — > K such that (fT3j) holds. For /i^ we get the following result: 
Proposition 3.5 pig fulfills the FKG inequality on R m ; i.e. 

[ H 1 H 2 d^> I Hidp^ [ H 2 dp h g V4ff 2 GO(l m ). (16) 

Clearly, (fT7)|) completes the proof of Theorem lH.il For Fi € 0(T), 1 — 1,2, one can 
choose the representing Hi from 0(M. m ) and from Prop. 13.51 and (|T5jl one obtains 
the correlation inequality for the measure defined on the right hand side of that 
equation for observables F t \, I = 1,2 and arbitrary A > 0, A C R d compact. 
Using Propositions (ii). (iii) and 13.21 (ii) one gets the convergence of both 

sides of the so-obtained inequality to the respective sides of (jEJ) with p = p g if we 
let first A / R d then A \ and finally carry out the short range approximation 
where in every step integrands in the numerators and denominators converge in 
L p (T,dpo). The statements Proposition I3.2ti3~4l (i) are not needed here, but in 
the following section. 

Proof of Prop. 13. 51 We first assume that the measure r defining ip is absolutely 
continuous w.r.t. to the Lebesgue measure ds on R with a Cq (R)-density function 
and suppr fl (0, oo) ^ 0, suppr H (— oo, 0) ^ 0. By investigation of the Fourier 
transform of p one gets: For / G C(R) exponentially bounded 

fdp = e-^ d J2 i -^T 1 I f dr * n 



n=l 



where r* n is the n-fold convolution of r with itself. 
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This implies that p is absolutely continuous w.r.t. the Lebesgue measure, 
dp(r)j) = g(r)j)dr)j with g E C 2 {R) and p > on 1, and log g is a well-defined 
function in C 2 (M.). Consequently, 



with W(rjj\j Eh) = U^(r]j\j Eh) — J2jeh diVj) ~ log^; — being the normal- 
ization constant J Rm e - ^'^ 61 ^ <8>jeL In order to verify the FKG inequality 
for a probability measure of type ()18jl with : R m — > R in C 2 (M m ), it suffices 
to verify the logarithmic FKG-criterium d 2 W(r)j\j Eh)/ drjjdrji < on IR m for 
I 3i lij £ L [ 1 9J . In our case 2nd order partial derivatives of W exist as 
v E C 2 (R) and log g E C 2 (M. d ). Carrying out the partial differentiations one gets 
from the definition of 



for j, I E L, j 7^ / and 0l defined in terms of the variables rjj as above. The 
crucial assumptions that v is concave, v" < 0, and G > enter only here. 

For the general case, where r is not necessarily absolutely continuous w.r.t. 
ds, we use the following approximation preserving the FKG inequality: Let dr e = 
[(1 - e)(xe * r) + e(x, * r e )] rfs where r e (A) = r(-A), A E B(R). x e C 2 (M) 
is a symmetric mollifyer, i.e. x > 0, x( — s ) — an d / R X^ — 1 ; Xe( s ) — 
x(es) / e. Then r e fulfills the above assumptions and furthermore 3C > such that 
supp r e C [— C, C] V0 < e < 1. Obviously r e — > r in law if e \ 0. Let p^ e be the 
lattice measure (fPfj) with r replaced with r e . To prove J Rm H dp^ t — > J Rm H dp^ 
as e \ it suffices to show the related statement for the free measures with 
v = as the multiplication of a continuous, exponentially bounded observable H 
with a Gibbs factor e~ u s again gives an observable of this type. As the single 
spin sites then are decoupled, we can identify L with {1, . . . , m} without loss of 
information. Thus, by (fTSj) . 



The integrals on the right hand side converge to the related integrals with e 
dropped as e \ 0: r e — > r in law implies r* n — > r* n in law and the support of r* n 
is contained in [—nC,nC] for 1 > e > 0. As if restricted to x 1 JL 1 [—NjC,NjC] is 
bounded and continuous, the convergence of the integrals follows. 

In order to show that the r.h.s. of (|2U|) converges to the left hand side of 
(fT9l with e dropped, it only remains to establish a summable upper bound. We 



notify that \H(s u ...,s m )\< Ke K ^ + - + ^ < Ke KC ^+- +N ^ for (s u . . . , s m ) E 



(18) 




(19) 




Weakly attractive interaction 



13 



supp 



< e < 1. Replacing the integrals on the r.h.s. of (J2Uj) with 

the latter e-independent term, we obviously get a finite sum. This finishes the 
proof. ■ 



4 The thermodynamic limit for the grand canon- 
ical measure 

Nelson's strategy for the proof of the existence of the thermodynamic limit of 
P{4>)2 Euclidean quantum field theories consists of two basic elements: Mono- 
tonicity and upper bounds for certain expectation values [TH]. The FKG in- 
equalities, established in the previous section, give us both elements of Nelson's 
strategy for the class of models under consideration, as we will prove below. First 
we however require one further technical assumption: 

Assumption 4.1 The linearly bounded, concave function v G C 2 (M) defining the 
interaction is monotonically falling, i.e. v' < 0. 

The point of Assumption 14.11 is that it renders — U g to be monotonically in- 
creasing in the following sense: From G > and (j3J) it follows that for rj, 7 G T, 
r) < 7 (i.e. r)(A) < 7(A) VA G B(R d )) G * rj < G * 7 -U g (j}) < -U B (i). 

The Assumption 14.11 at the first look might look more restrictive than it is: 
Given a v G C 2 (R) linearly bounded and concave that is not monotonically falling, 
one can do the following replacements 

v(4>) -> v((f))-b(f), &>|H|oc, 

e -sb(3\\G\\ L i (M d idx) 

dr<3) " J^^m^^ dr(s) ' (21) 

z -> z [ e - sb ^ G hHmi, d *) dr(s) . 
J[-c,c\ 

Heuristically, the re-defined system, which now fulfills Assumption I4.1[ and the 
original one describe the same physics. The only difference is that a linear term, 
i.e. a self-energy term proportional to the charge s of the particle, see ((H), has been 
subtracted from the interaction potential and has been added to the chemical po- 
tential. In our case, where we deal with particles of variable charge s G [—C,C], 
this means that on the one hand we have to re-define the relative chemical po- 
tentials of particles with different charges, i.e. to re-define r, and to re-define the 
over-all activity z on the other hand. This explains equation (|21j). 

A more technical issue however is touched by (f2Tj) which in a large sense can 
be seen as the issue of boundary conditions. To explain this, let us for simplicity 
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chose r = (l/2)(5+i + 6-1) and that v is concave with v (—(/>) = v(4>). The system 
then is formally invariant under the replacement 77 — > — 77. The re-definitions 
(|2ip break this invariance if the infra-red limit g S (3 has not been taken yet: 
The effect of the re-definition of the energy density v leads to a gain of potential 
energy of the "+" particles w.r.t. the "— " particles only in the neighborhood of 
supp g, g being the IR-cut-off function, whereas the gain in the relative chemical 
potential for "— " charged particles takes place everywhere. Hence, outside the 
support of g, negatively charged particles dominate. This can be interpreted as 
boundary conditions of "—"-type. 

A related discussion with Assumption 14.11 formulated for v monotonically in- 
creasing and b in (|21j) replaced with —b would then lead to "+"-type boundary 
conditions. For data leading to a phase transition (in the language of the fer- 
romagnetic spin systems of Section 3 this would be spontaneous magnetization) 
"+" and "— " type boundary conditions probably cause different phases. But a 
detailed discussion of this issue is beyond the scope of this article. 

From now on we consider only energy densities v such that Assumption 14.11 
holds. At first we want to establish monotonicity of expectation values in g for 
monotonically increasing function: 



Proposition 4.2 Let F e 0(T) and g 1 < g 2 , g u g 2 G C (R d ), g u g 2 > 0. Then, 



F dfi gi < Fdfi 



'9-2 



(22) 



Proof. It suffices to show that for F G 0(T), g,f e C (R d ), g,f>0, we get 
(d+ Jr F dfi g /df) = lim <Xo (/ r F d^ g+tf - J r F d/i g )/t > 0. 

We note that lim^ (e _f/9+t/ — e~ Ua )/t = —Ufe~ Ug pointwisely on T and that 
the differential quotient on the left hand side for < t < 1 has a L P (T, /i )-upper 
bound, e.g. b(\i]\, f + gjd'tMJ+g) _ Hence one gets by application of the quotient 
rule that 



d+ jr F dn 9 
df 



0, 



J T Fe Ua dfJL 



dfio 



df 



F(-Uf)dn g - Fd^i 



-U f ) d\i 



9 ■ 



(23) 



We notify that by Assumpt ion 14 . 1 1 ( — Uf) formally is of monotonic increase, hence 
the right hand side of ()23|) formally is non-negative by Theorem 13.11 as required. 



To make this argument rigorous, one could approximate {—Uf) with observ- 
ables in 0(T). Here we use the lattice approximation instead. First we recall 
that (-Uf) e 0{R m ) by Assumption O Hence, the right hand side of (J23J) 
is non-negative when going to the lattice, cf. Proposition 13.51 By Propositions 
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I3.2H3.4I products of terms like e u s,& } (—Uf A ) and F A converge in L p (T,du, ) to 



the corresponding products of e Ua , (—Uf) and F when first A /* M. d , then A \ 



f) 

and finally an approximation of G by functions of Co(M d ) is used. Thus, the right 
hand side of (J23"j) can be approximated with non-negative expressions and is thus 
non-negative. ■ 

The next proposition gives upper bounds by comparison of the interacting 
ensemble with a gas of non-interacting particles with a different (space-dependent) 
activity. Let yUo, 9 be the measure on (T,B(T)) obtained by replacing v(<p) in the 
definition of ri g with the linear function —b(f). 

Proposition 4.3 Let F G 6(T). Then, 

J Fdu. g < j Fdfio, g <M VgeC (R d ),0<g<(3 (24) 

where M is finite and depends only on F , z, C, b, ||G|| L i( K d jrf:r ) and (3. 

Proof. For a G [0, 1] let fi a>g be the measure obtained by replacing v(4>) in the 
definition of \i g by av((p) — (1 — a)b<p. This density is again concave for a G [0, 1], 
hence \i a%g fulfill the FKG correlation inequalities. Obviously fi a<g interpolates 
between /x ,q an d ^i, g = n g - In order to prove the first inequality in ([23$ it is 
sufficient to prove d LF du. ag /da < 0. 

Let U a! g be the potential corresponding to the energy density av((p) — (1 — a)b(j) 
and let U' = dU atg /da, i.e. U' is the potential with IR-cut-off g correspond- 
ing to the energy density v(4>) + b<p. By the same reasoning as in the proof of 
Proposition 14. 21 one gets 



dj T F dri at: 
da 



J F(-U' a>g )d» a , g - J F d^ g J(-U' a>g ) d^ g . (25) 

We note that U' a g is of monotonic increase in the same sense as explained for 
—Uf in the proof of Proposition 14.21 By the lattice approximation Prop. I3.2H3.4I 
and the FKG inequality on the lattice Prop. 13.51 one thus gets as in the previous 
proposition that the r.h.s. of (|23|) is non-positive. This proves the first inequality 
in (jUJ). 

To establish the uniform bound M, we note that one gets by a standard 
calculation on Laplace transforms (see e.g. 0) 

(MMHVJ) dfiQ ^ = // E 4 I - C , C ][<= M " W+s/,! ' L l]*(^!/ ) (26) 



r 



/, h G L 1 (IR d ,(ix) bounded and continuous. We set / = bG * g (i.e. —U 0tg (rj) 
(r], /)) and we chose h G C f . f .(K d ), h>0, such that F(rj) < Ke {M ' h) . Then, 

J re (\r,\,h)+(v,bG*g) dfiQ ^ 

Fd ^°> 9 ~ K j r e^ bG ^du. ( v ) 



r 



16 



H. Gottschalk 



= ^/ Rd / [ _c,c][ e!sl ' l(! ' )+s,,G * 9(! ' ) - eS,,G * 9(! ' ) H'-W^ (27) 

For < g < ft, \\G * g\\oo < /3\\G\\ v-md^y The integrand in the exponent on the 

r.h.s. thus is smaller than Rh(y) G L^R^cix) with R = Ce CMoo+cmmLl ^ d > d *) . 
■ 

We are now in the position to formulate the main theorem of this article. The 
notation g /* ft which is needed to formulate the infinite volume limit has the 
following precise meaning: For any sequence {g n }n G N of non- negative functions 
in C (lR d ) such that for any N G N 3n G N such that g n (x) = ft for all x, 
\x\ < N, and n > no the limit exists and is unique. 

Theorem 4.4 Let g / ft. Then, 

(i) There exits a uniquely determined measure fip on (S'(M. d ),B) such that fi g — > 
fi/3 weakly. 

(ii) up is invariant under translations. If G is invariant under rotations, fip is 
also rotation invariant. 

Proof, (i) We notify that by Propositions I4.2H4.3I lim g yp j r F dfi g exists for 
F G 6{T). 

Let now F G 0(Y) such that F(rj) = H({r), h x ), (77, h 2 )), h u h 2 G C f . f .(M d ) and 
H : M. 2 — > R in C 2 (M. d ) together with first its and second order partial derivatives 
exponentially bounded. Then \im g yp J r Fdfi g still exists. In fact, one can write 
F = F^ — F^- with F^ G 0(T). As the limit exists for F^ separately, it also 
exists for F. One choice of the F*^ is as follows: Let H ± be the positive/negative 
part of a function H = H + — H . We set 

H u (x,y) = #±(0,0)+ [\d 2 H(0,t)) ± dt 



Jo 

(d 1 H(s,0)) ± ds+ / / (d l d 2 H(s,t)) ± dsdt (28) 
Jo Jo 

for x,y G 1R, di/ 2 the partial derivative w.r.t. the first/second argument and 
Fn(r l ) = Hn({ Vj h 1 ),{r ] ,h 2 )). 

In the next step we prove weak convergence to a measure on (<S'(M d ), £>)) using 
Minlos' theorem |H]: Let / G S(R d ) then /+, /" are non-negative functions in 
C { ,(R d ). Then, F^) = cos(( V , f+) - (rj, /")) and F 2 ( V ) = sin(( V , /+) - ( V , /")) 
are observables of the type described in the preceding paragraph. Hence the limit 
of the characteristic functional C g (f) = J r e l( ~ v '^dix g (r]) converges to Cp(f) G C 
for / G S(R d ) as g / ft. 

Cp : 5(lR d ) — >• C is positive definite and normalized as the limit functional of 
functionals with these properties. It remains to prove that Cp is continuous w.r.t. 



Weakly attractive interaction 



17 



the S (R d ) -topology. For fixed cut-off g one gets 

\c g (f)-c 9 (h)\ < v [(v,(f-hyyd» 9 (v) 



< 



£ 

<T=± 



(77, (/ - hf)- dfi (v) + / (V, (f ~ hY) + d/jLoM 
J Jr 



(29) 

In the last step we applied Proposition 14.21 and 14.31 to the decreasing observables 
(77, (/ — h)^)' and the increasing observables (77, (/ - /i) ± ) + . As h -> / in S(M d ) 
it follows (f — h) ^ — > in Cf.f(lR d ) and the rf/io-integrals vanish under this limit, 
as can be seen by dominated convergence in L p (r, //o)- 

We have to give a g- independent estimate for the vanishing of the /i 0i9 -integrals 
for < g < (3. One can proceed as follows: 

(77, (/ - h)±) + dptoM < [ (e^ f - h ^ - 1) dfioM 
v Jr 



as (h — /) + — >• in Cf.f(R ) with R = R((h — f) + ) as in the proof of Proposition 
14.31 independent of g and bounded for (h — — > in Cf,{{R d ). 

We have proven that Cp is a characteristic functional and \ip can now be 
defined as the unique measure on (<S'(R d ),23) with the given characteristic func- 
tional, p3]. As weak convergence is equivalent with the convergence of charac- 
teristic functionals by Levy's theorem, the assertion of (i) follows. 

(ii) Let G be invariant under rotations and reflections. For rj E T, 77 = 
Syesuppr; s y$yi x E R d and D an element of the orthogonal group on M d let 
V{x,d} = J2ye S up P r, s y S Dy+x and KM = h(D~ 1 (y - x)), h : R d -> R. For 
F E L 1 (r, / u ) let T{ XjD }F(i]) = F(r]{ Xt Dy). /i is invariant in distribution under 
Euclidean transformations, i.e. Tf d \^jlq = fMo with T? x D i the dual action of T^.z?} 
(this e.g. can be deduced from the invariance of (J2J) under such transformations 
and the uniqueness statement of Minlos' theorem). For the potential energy 
we get U g {r] {x , D} ) = j Rd v(G * n {x>D} ) gdx = j Rd v{{G * r]) {x , D} )gdx = j Rd v(G * 
V) 9{x,D}- 1 dx. Thus, Tf D ^ji g = ^ 9{xDy The assertion now follows from the 
equivalence of the limit g /* [3 and g{ x ,D] S ft- If G is not invariant under 
rotations, the argument still holds for D — 1. ■ 
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